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Abstract 


The  Jaumann  stress-strain  approach  has  been  used  to  evaluate  a 
nonlinear  structural  response  to  a  shell  like  geometry.  The  method  employs 
a  finite  element  solution  incorporating  the  Jaumann  stress-strain 
relationships  based  on  large  displacement  and  large  rotation  using  a 
corotational  technique.  The  resulting  equations  include  the  continuity  of 
stresses  and  displacements  in  the  thickness  direction.  Thus,  it  is  possible  to 
evaluate  the  effects  of  direct  and  shear  stresses  within  a  laminated 
structure. 

Initial  comparisons  of  the  method  of  analysis  has  been  made  with  two  other 
theories.  The  first  is  the  total  Lagrangian  theory  based  on  an  in-house 
program,  and  the  second  is  an  Eulerian  based  theory  defined  by  ABAQUS. 
Several  problems  have  been  attempted  including  a  cylindrical  shell  panel 
made  from  composite  materials  undergoing  nonlinear  collapse.  Two  new 
algorithm  have  been  developed;  the  first  to  consider  a  nonconservative  force 
system  and  the  second  to  transform  the  resulting  equations  from  a 
curvilinear  coordinate  system  into  a  Cartesian  system.  The  last  algorithm 
defines  a  global  set  of  equations  which  can  be  used  for  the  mixture  of  various 
element  geometries. 

Results  show  that  the  need  for  a  correct  through  the  thickness 
evaluation  of  the  stresses  of  the  shell  is  required  as  the  shell  thickness 
increases  relative  to  the  resulting  geometry.  In  a  predominantly  membrane 
resistant  structure,  the  three  types  of  theories,  Eulerian,  Jaumann  and 
Lagrangian  show  close  comparison. 
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1  Introduction 


The  focus  of  this  research  is  to  carry  out  an  analysis  of  satellite  antennas 
made  from  various  composite  materials  incorporating  a  program  which  can 
be  specifically  directed  toward  this  type  of  structure  without  resorting  to 
commercial  programs.  In  this  inflatable  structure,  a  paraboloidal  shape  is 
desired  for  optical  and  RF(Radio  frequency)  reflectors  [l][4][5]  because  of 
their  scientific  use.  Thus,  the  major  concern  is  to  obtain  a  membrane  space 
structure  which  approaches  a  paraboloid  geometry  under  internal  pressure. 
Therefore,  this  research  is  an  attempt  at  developing  a  robust  but  efficient 
finite  element  program  which  can  simulate  the  major  characteristics  of  the 
need  for  arriving  at  an  approximate  final  geometry.  The  method  employed  is 
one  which  incorporates  the  concept  of  large  movement  due  to  a  nonconserva¬ 
tive  pressure  loading  (one  that  remains  normal  to  the  deforming  structure). 
This  method  is  based  upon  the  Jaumann  stress-strain  relationships  using  the 
corotational[2]  [3]  technique.  The  resulting  equations  include  the  continuity 
of  stresses  and  displacements  in  the  thickness  direction.  Thus,  it  is  possi¬ 
ble  to  evaluate  the  effects  of  direct  and  shear  stresses  within  a  laminated 
structure  and  trace  any  nonlinear  geometric  characteristics. 

In  this  first  report,  the  Jaumann  relations  are  developed  within  a  total 
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Lagrangian  system  using  a  corotational  approach.  An  in-house  computer 
program(JAGS1)  is  modified  to  handle  the  appropriate  geometries  and  results 
are  compare  to  an  Eulerian  based  program,  ABAQUS.  Several  problems  have 
been  attempted  including;  a  cylindrical  shell  panel  made  from  composite 
materials  undergoing  nonlinear  collapse;  an  isotropic  membrane  acted  upon 
by  a  nonconservative  force  system  and  a  parabolic  shell  acted  upon  by  a 
pressure  loading.  These  problems  verify  the  accuracy  of  the  theory  and 
program. 

Results  show  that  the  need  for  a  correct  through  the  thickness  evaluation 
of  the  stresses  of  the  shell  is  required  as  the  shell  thickness  increases  relative 
to  the  resulting  geometry.  In  predominantly  membrane  resistant  structures 
the  comparison  among  the  three  types  of  theories,  Eulerian,  Jaumann  and 
Lagrangian(SLR2  [6] ,  LDLR3[7])  are  very  close. 

2  Theory 

In  this  section,  a  brief  review  of  the  corotation  theory  leading  to  the  Jaumann 

stress-strain  relationship  is  carry  out.  This  theory  is  presented  to  show  a  total 

1  Jaumann  Analysis  of  General  Shells 

2  Simplified  Large  Displacement  and  Moderately  Large  Rotation 

3  Large  Displacement  and  Large  Rotation 
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Lagrangian  corotational  approach.  Much  of  this  development  can  be  found 
in  several  papers  written  by  Palazotto  et.  al.[7]  -  [10], 

2.1  A  Total  Lagrangian  Corotational  Finite  Element 
Scheme  used  in  the  Analysis  of  An  Inflatable  Shell 


Figure  1:  Deformation  of  an  infinitesimal  volumn  element. 
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The  theory  makes  use  of  the  polar  decomposition  method  to  facilitate 
the  use  of  a  local  (and  linear)  displacement  field  at  an  infinitesimal  region  of 
interest  in  the  nonlinear  deforming  body.  The  Jaumann  stress  [11]  Jmn  and 
strain  Bmn  are  given  by 

Jmn  =  2  dxx  dx2  dx3  (dXm  fm  ' ln  +  dXn  fm  ' lm)  (X) 


Bmn  - 


1  /  du 


2  \dxm  "  dxn 


where  the  fm  are  the  force  resultants  on  the  faces  of  the  deformed  paral¬ 
lelepiped  (  Figure  1).  For  example,  fi  acts  on  the  deformed  dx2  -  dx3  plane. 
The  local  displacement  vector  u  of  equation(2)  is  measured  with  respect  to 
the  displaced  location  of  a  material  point,  hence  at  any  given  point  on  the 
deformed  reference  surface,  u  =  0,  though  its  derivatives  (which  will  give  rise 


to  the  strains)  are  non-zero.  The  Jaumann  stresses  and  strains  are  defined 


with  respect  to  the  orthogonal  directions,  denoted  by  unit  vectors  ifci  associ¬ 
ated  with  the  stretched  and  rigidly  rotated  volume  element  (Figure  1(a)).  On 
the  other  hand,  the  second  Piola-Kirchhoff  stresses  Smn  and  Green’s  strain 
Lmn  are  defined  by 


dx^  dx 2  dx3 


dx^m)  f(m)  =  ( Sm[n )  i(n)  j 
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where  An  is  the  magnitude  of  so-called  “lattice  vector”,  and  the  parentheses 
suspend  the  tensor  summing  convention.  In  the  figure,  it  is  seen  that  the 
Second  -  Piola/Green  measures  are  associated  with  the  directions  along  the 
deformed  (and,  in  general,  not  orthogonal)  edges  of  the  element,  as  shown 
in  Figure  1(b).  The  directions  of  the  lattice  vector  correspond  to  the  direc¬ 
tions  of  the  unit  vectors  i^.  So,  in  general,  the  component  of  the  Second 
Piola  stresses  are  along  neither  those  of  the  undeformed  coordinate  system, 
as  are  the  Eulerian  measures  depicted  in  Figure  1(c),  nor  those  of  its  rigidly 
translated  and  rotated  counterpart  in  the  deformed  body  (as  are  the  Jau- 
mann  measures).  This  is  a  consequence  of  the  Green  strains  being  energy 
related  measures  rather  than  strictly  geometric  measures,  like  Jaumann  or 
engineering  strains.  To  use  the  Jaumann  measures,  which  are  local,  the  effect 
of  rigid  body  translation  and  rotation  must  be  removed  so  that  the  effect  of 
stretching  is  seen. 

In  a  layered  composite  consisting  of  N  layers,  the  local  displacement 
vector  with  respect  to  the  local  £ri(  coordinate  system  of  Fig.  (2)  as  presented 
in  Pai  and  Palazotto[3]  is  defined  as 

u  =  u(il)ii  +  u%h  2  +  t4l)i3  (5) 
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Figure  2:  Infinitesimal  element  undergoing  deformation(after  Pai  and  Pala- 
zotto,  1995) 
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where 


ui]  =  ui(x’y)  +  ^[^(x,y)-9°2(x,y)}+l5z  +  a<t)(x,y)z2+P?(x,y)z3 
u*  =  U2(x,y)  +  z[e1(x,y)-e°1(x,y)}+'r 4z  +  a^(x,y)z2  +  ^(x,y)z3 

u3}  =  «3  (^, 2/)  +  4°  (x,  y)z  +  0®  (x,  y)  z2 

Here,  w°  (j  —  1,2,3)  are  the  components  of  displacement  (with  respect  to 
the  local  coordinate  system  £r?C)  of  a  point  which  is  located  on  the  reference 
surface  at  (x,  y)  before  deformation.  The  rigid  body  rotations  and  shear  ro¬ 
tations  are  given  by  6  and  7,  respectively.  The  angle  between  the  transverse 
coordinate  2  and  the  normal  to  the  reference  surface  in  the  undeformed  con¬ 
figuration  as  measured  in  the  rz  plane  is  given  by  0°.  The  corresponding 
angle  in  the  yz  plane  is  given  by  6^-  The  shear  rotation  angle  in  the  xz  plane 
at  the  reference  surface  is  denoted  by  7S,  and  represents  the  rotation  of  the 
normal  to  the  reference  surface  due  to  transverse  shear  deformation.  The  cor¬ 
responding  angle  in  the  yz  plane  is  y4.  The  terms  0$  and  /?«,  are  referred 
to  as  shear  warping  and  thickness  stretch  functions.  These  functions  are 
used  to  describe  the  kinematic  behavior,  beyond  simple  rotation  of  the  rigid 
normal,  of  the  material  away  from  the  reference  surface,  and  allow  coupling 
of  the  displacement  u[ )  and  1$  via  the  shear  angles  at  the  reference  surface. 
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That  is,  75  can  affect  displacement  through  the  warping  functions  and, 
likewise,  74  can  affect  u%\  By  defining  the  shear  warping  functions,  Gi  and 
G2,  and  the  thickness  stretch  function,  G3  as 


=  75  z  +  ai^+Z^z3,  G2  —  74z+c4,)22-l-/?;;,)23and  G3  =  a^z+fflz2, 


(6) 


the  kinematics  of  the  eqn.(6)  may  be  written  as 


ui  ui  (xi  y)  +  z  [$2  (£>  y ) —  $2  (®i  2/)]  ■+•  Gi 
u2]  =  ^2  (^, y)  +  ^  [<9i  (a;, y)  —  <9? (x, y)]  +  G2 
uz]  =  ul(x,y)  +  G3.  (7) 

The  stretch  function  G3  is  usually  small,  especially  for  thin  shells  and  one 
may  neglect  G3  and  its  derivatives  in  most  strain-displacement  expressions. 
This  is  based  upon  the  claim  that  the  effect  of  transverse  shear  strain  on 
the  in-plane  strain  is  negligible.  Under  this  assumption,  the  Jaumann  strain 
displacement  relations  become 


du  . 

p>\z) 

dx  11 

du  . 

b22 

=  TP  •  12 
dy 

dW 

-°33 

du 

=  dl‘13 

^  •  ii  =  (1  +  ex)  cos  761  —  1  +  2  ( \  -  fc°)  +  Gitx  -  k5G2  (8) 
•  i2  =  (1  +  e2)  cos  762  -  1  +  2  (k2  -  k%)  +  G2<y  -  kAGx  (9) 
*3  =  1  (10) 
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(11) 


2  Bw 


23  — 


2  Bg 


2B(i) 


12  — 


du  .  du  . 

Qy  ’  *3  +  '  12  =  W2]Z  —  k$2G\  —  ^2^2 

du  .  5u 

dx  13  **"  ~dz  =  ^1,z  —  —  ^61^2 

du  .  du  , 

a^',2  +  ^'“ 


-  (1  +  ei)  sin  761  +  (1  +  e2)  sin762  +  2:  (fc6  -  kfj 

+C?2,z  ~  ^62^1  —  k2G2 


(12) 


(13) 


where  k6  =  k61  +  k62  and  k%  =  A&  +  k%2. 

For  comparison,  the  total  Lagrangian  LDLR  making  use  of  the  Green 
stain  tensor,  incorporates  the  kinetics,  written  in  global  Lagrangian  cylindri¬ 
cal  coordinates,  as 


«i  =  tt  +  ^sin(^1)cos('02)  (14) 

u2  =  v(l-^j+zsin(ip2) 
uz  =  tw  +  ^{cos(^1)sin(^2)  —  1} 

where  ipx  is  the  angle  made  by  the  deformed  normal  to  the  original  normal 
in  the  x-z  plane,  and  ip2  Is  the  angle  made  by  the  deformed  normal  to  the 
original  normal  in  the  y-z  plane.  The  bending  angles  ipx  and  ip2  are  similar 
to  Euler  angles  for  rigid  body  motion,  u,  v  and  w  are  the  displacements  of 
i-th  layer  on  the  mid  plane  and  R  is  the  radius  of  the  shell.  This  is  a  first  or¬ 
der  shear  theory  with  a  linear  through  the  thickness  shear,  and  the  kinematic 
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assumption  violates  the  zero  shear  condition  on  the  top  and  bottom  surfaces. 
Also,  it  can  be  observed  that  by  using  a  small  angle  approximation  for  both 
ip i  and  •02)  the  basic  kinematics  used  in  the  SLR  theory  is  obtained.  The 
kinematics  used  in  the  SLR  theory  are  based  on  the  third  order  shear  defor¬ 
mation  theory  which  assumes  that  the  through  thickness  shear  is  a  parabolic 
curve.  This  theory  also  makes  use  of  the  Green  strain  tensor  but  assumes 
relative  small  rotations.  A  detailed  description  regarding  the  derivation  of 
the  applicable  kinematics  that  can  be  used  for  large  kinematics  is  presented 
in  Gummadi  and  Palazotto[13]. 

The  relationship  between  the  Jaumann  stresses,  Jmn,  and  the  Jaumann 
(or  Biot-Cauchy  Jaumann)  strains,  Bmn,  for  the  ith  lamina  of  transversely 
isotropic  material  may  then  be  written  as  (for  the  transformed  relationships) 

51?  52  52  52  o  oil  b<? 

Ql2  Q22  ^23  Qm  0  0  B'22 

52  52  52  52  0  0  j?g 

1  1  •  (15) 

52  52  52  52  0  0  2  Bg 

0  0  0  0  52  5*2  2  Bg 

0  0  0  0  ZJ2  52  2Bg 
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The  first  variation  of  potential  energy  on  an  element  basis  is  represented  as 

™  =  5g//jC(*[B«]r  Qw]  [B«]  +  [B'f  *J®]«[B«])-V 

=  J  J  J  $  [B(t)]  dV  (N  =  number  of  layers)  (16) 

where  V  is  the  undeformed  volume  of  the  shell  element  structure  and  (i) 
refers  to  the  value  of  the  function  in  the  ith  layer  of  the  laminate.  And  the 
Q(t)  is  a  transformed  stiffness  matrix  of  principal  stiffness  matrix  for 
the  ith  laminar  with  respect  to  the  £  axis,  while  B(<)  is  the  Jaumann  strain 
in  the  ith  laminar  with  respect  to  the  £  axis.  Figure  (3)  shows  a  four  noded 
44  degree  of  freedom  (DOF)  finite  element.  The  11  degree  of  freedom  at  each 
corner  are  u ;  u,x ;  u,y ;  v;  v,x ;  v,y ;  to;  w,x ;  w,y ;  74  and  75.  Hermitian  shape 
functions  are  used  for  all  DOF  except  the  transverse  shear  DOF  74  and  75, 
which  use  bilinear  shape  functions.  In  order  to  solve  the  nonlinear  finite 
element  equations,  the  Newton-Raphson  method  is  used 

|  M  W*}  - 1  (W  -  M  R})^,  an 

where  |R^j  is  the  elemental  nodal  loading  vector,  Jk^'I]  is  so  called  tangent 
stiffness  matrix,  {q0}  is  the  displacement  vector  at  the  last  converged  incre¬ 
ment.  jqblj  represents  the  current  nodal  displacements,  and  Ne  is  the  total 
number  of  elements  making  up  the  model.  The  product  [k^]  {qbl}^  ^  ^ 
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is  a  vector  described  as  follows, 


M  W,V,H<„  -  SSm  ifjjf ) 


(18) 


The  formulation  of  the  matrices  in  the  above  equation  is  described  in  detail 
in  the  reference  of  Pai  and  Palazotto[ll]  also  in  Greer  and  Palazotto[12]. 
Detailed  explanation  for  the  terms  in  equations  (18)  are  shown  in  Appendix 
A.  These  terms  relate  to  the  variation  of  the  functions  of  the  displacement 
gradient.  In  addition  ['P°]T'  is  a  known  matrix  at  fq0} 

24x12  1  J 


2.2  Global  Coordinates  Representation  of  Elements  Us¬ 
ing  Transformation  Matrix 

In  order  to  incorporate  an  element  developed  in  a  local  curvilinear  coordi¬ 
nate  system  with  other  elements  in  their  own  local  coordinate  system,  one 
must  carry  out  a  transformation  to  a  global  coordinate.  It  is  appropriate 
to  make  use  of  a  Cartesian  system  for  this  purpose.  In  order  for  this  to  be 
establish,  one  must  evaluate  a  transformation  based  upon  the  appropriate 
orthonormal  system  and  its  direction  with  respect  to  the  Cartesian  axes.  In 
this  section,  the  authors  develope  the  relationships  between  the  Jaumann 
coordinate  system,  employing  a  local  curvilinear  coordinate  system  and  the 
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global  Cartesian  coordinate  system.  (As  yet,  these  expressions  have  not  been 
implemented  within  the  program) 

It  is  possible  to  express  the  equilibrium  state  in  the  original  local  coordi¬ 
nate  system  using  the  equation 


[Ka\  K}  =  {i?Q} 


(19) 


where  subscript  a  represents  local  curvilinear  coordinate  systems. 

If  we  were  to  express  the  same  state  of  equilibrium  in  a  Cartesian  system, 
one  would  find 


m  {ux} — {Rx} 


(20) 


where  subscript  x  represents  global  Cartesian  coordinate  systems. 


The  transformation  matrix  [T]  between  the  local  curvilinear  coordinate 


system  and  global  Cartesian  coordinate  system  for  an  element  is  next  intro¬ 
duced  yielding 


K}=  [■ T ]  {«*} 

44x1  44x1 


(21) 


The  [T]  matrix  is  a  matrix  that  has  a  rank  of  44x44  if  one  considers  a  shell 
element  as  in  Figure  (3).  Thus  each  degree  of  freedom  in  the  Cartesian  system 
is  transformed  into  the  appropriate  degree  of  freedom  in  the  curvilinear 
system.  The  [T]  matrix  is  developed  as  following;  The  vector  A  can  be 
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Figure  3:  44  DOF  Shell  Element 

expressed  using  direction  cosines.  In  Figure  (4),  the  vector  A  is  drawn  from  an 
arbitrary  point  P  in  space  and  will  be  referred  to  a  set  of  unprimed  Cartesian 
coordinate  axes  with  the  origin  at  O,  and  to  a  set  of  primed  coordinated  axes 
with  the  origin  at  O'.  Both  sets  of  axes  may  be  translated  without  rotation 
to  the  common  origin  at  P.  It  is  more  convenient  to  summarize  these  sets  of 
direction  cosines  in  tabular  form  as  shown  in  Table  1. 


19 


Figure  4:  Representation  of  a  vector  in  two  sets  of  coordinate  axes  with  a 
different  orientation 
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Table  1:  Direction  Cosines  between  two  coordinate  axes. 
Then,  the  transformation  matrix 
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m3 

(22) 
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yielding 
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Now,  the  nodal  degree  of  freedoms  of  one  of  the  four  noded  element  is 
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Because  of  the  fact; 


du  _  dudx  du  dy  du  dz 

dx'  ~  +  +  lU’y+n'u’*  (25) 

h  mi  0 

du1  _  du'  dx  du'  *dy  'duTdz' 

dx 1  dx  dx1  +  dy  dx'  +  ~dzdx'  (26) 

=  l\{l\U  +  miv  +  niw),x  +mi(liu  +  miv  +  n\w),y 

=  hhux  +  rriiliu,y  +m1liv,x  +m1m1u,J/  +nll1w,x  +n1m1w,3, 
dv  _  dv  dx  dv  dy  dv  dz 

dx'  dx  dx'  +  dydx'  +  ~dz~dx'  =  ll  V,x  +mi  v,y  +niV’z  (27) 

etc. 

The  transformation  matrix  can  be  expressed  for  the  44  DOF  element  as 
follows 

S'  [0]  [0]  [0] 

[0]  S'  [0]  [0] 

[0]  [0]  S'  [0] 

[0]  [0]  [0]  S' 

where  [0]  is  11  x  11  null  matrix  in  which  all  the  component  of  matrix  are 
zero.  Equation  (19)  can  now  be  expressed  as 


[jq  [t\  K}  =  {*„} . 


(29) 


Multiply  both  sides  of  the  equation  (29)  by  [T]  1  =  [T]t  yields 


where 


mw  =pf  {*»}  • 


"V- 

Ka 


V“ 

Rx 


(30) 


K,  =  [7f  (*•„]  [T]  and  ii,  =  [Tf  {/(„(. 

The  Jaumann  strain  {B}  is  a  function  of  ua  which  is  {£}  =  {B  (ua)}.  It 
can  be  expressed  in  global  Cartesian  coordinate  system  using  relationships 


B'ij  =  UkljiBu 

(  hj  —  direction  cosines  between  the  local  and  global  coordinate  systems). 

And  Jaumann  stress  J' , which  is  function  of  Jaumann  strain  B'  can  also 
be  expressed  as 


J'ij  —  hkljiJki- 


2.3  Nonconservative  Loading  Case  in  JAGS 

It  is  required  that  the  loading  system  considered  in  an  antenna  analysis  be 
directed  normal  to  the  deformed  surface.  Thus,  for  this  to  occur  one  must 
consider  nonconservative  loading  equations  for  evaluation.  The  loading  must 
be  applied  in  increments  normal  to  the  surface. 
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First  express  the  stiffness  matrix  for  the  ith  load  increment  as; 

(Jf„)  =  [iC(]  +  (31) 

note  [Ki]  can  be  a  nonlinear  stiffness  matrix,  and  [A K\n  =  dKn  =  . 

(For  example  Kn  =  3  (u?  +  u%)  then  dKn  =  3  (2u1du1  +  2u2du2).  The  equi¬ 
librium  equation 

[K]  M  =  {R}  (32) 

can  be  replaced  by 

([K]  +  [Aif]J  (M  +  {Aun})  =  {R}  +  {A  (33) 

resulting  in 


[K]  {u}  +  [K]  { Aun}  +  [AKn]  {n}  +  [Aifn  ]  {Aun}  =  {i?}  +  {A R,}  (34) 

The  term  [Aif„]  {Aun}is  H.O.T.,  so  we  can  ignore  it.  Equation  (34)  reduces 
to 

[K]  {u}  +  [K]  {Aun}  +  [A Kn)  {u}  =  {Ai^}  +  {/?} .  (35) 

Since  [if]  {u}  =  {i?},  one  can  subtract  these  terms  out  of  the  solution  in  order 
to  find  a  resulting  new  incremental  equation  of  state.  The  Newton-Ralphson 
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Technique  can  be  used  to  solve  this  nonlinear  set  of  equations  giving 


[Kt  }  8  {An}  =  {A R*}  -  [K]  {A un}  +  [AKn]  {u}  (36) 

where  Kj  =  ([A/fn]  {u}  +  [K]  {Aun})  <5  {Aun}  (referred  to  as  the  tan¬ 
gent  stiffness  matrix).  The  terms  [Aif„]  {u})  are  evaluated  by  con¬ 

sidering 

^  =  [AK]{u}  +  [K}{Aun}-{ARn}=0 ,  (37) 

which  is  the  incremental  equilibrium  relation. 

There  are  i-equations  with  i-DOF,  then  the  expression  using  the  index 
notation  is 

ipi  =  AKij  Uj  and  d'ipi  =  K{j  uj  (38) 

where 

Kij  =  (A Kinlln)  ,j  .  (39) 

Therefore 

Kt  =  [( AKin  un)  ,j  +K]  8  { Aun}  (40) 

For  example,  let’s  consider  the  simple  two  dimensional  case, 

K11U1  +  K12U2  =  R\; thus  ip1  =  K\\ u\  +  K12U2  —  Ri  —  0  (41) 
K21  ux  +  K22  u2  =  i?2 ;  thus  ip2  ~  K21  ux  +  K22  u2  -  R2  =  0  (42) 
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Then, 


fo/h  _  g(gn«i)  ,  9 {K12 U2)  _  d(Kinun)  d^i  _  d (ifin un) 
du\  du\  du\  du\  ’  du2  du2 

djh  =  d(K2i  t*Q  ^ (^22 _  d (K2n un)  dj>2  _  d (K2n un) 
du\  du\  du\  du\  ’  du2  du2 

If  we  use  the  modified  Newton-Raphson  technique,  then  K  6 A uni+l  = 
A  Rn  -  KAuni  +  AKniUi  can  be  used  to  evaluate  6Au„i+1.  The  value 
of  Au„i+1  (normal  increment  of  displacement  for  the  i  +  1  increment)  = 
Auni+6Auni+l.  The  satisfaction  of  the  incremental  state  of  equilibrium  is 
carried  out  within  an  established  tolerance. 


3  Application 

Several  attempts  have  been  made  to  compare  the  method  of  the  analysis  of 
the  total  Lagrangian  theory  based  on  JAGS,  and  the  Eulerian  based  theory 
defined  by  ABAQUS.  Three  cases  have  been  performed.  The  first  case  was 
cylindrical  shell  acting  under  a  concentrated  transverse  load  as  in  Figure  5. 
This  case  depicts  a  shell  panel  undergoing  movement  which  includes  an  insta¬ 
bility  or  collapse  point.  Thus,  geometric  nonlinear  displacement  and  rotation 
can  occur.  This  usually  is  displayed  in  a  problem  of  this  form  by  showing  a 
load  vs.  displacement  curve  for  the  degree  of  freedom  at  the  applied  force. 
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There  are  two  regions  for  this  type  of  curve.  The  rising  or  stable  region 
and  the  descending  or  unstable  region.  Each  point  in  the  plot  represents  an 
equilibrium  state. 


Figure  5:  Simply  supported  point  loaded  cylinderical  shell 


The  material  properties  are  described  as  in  Table  2.  The  shell  properties 
were  modeled  using  a  Hercules  AS4-3501-6  graphite  epoxy  composite  with 
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Ex 

E2 

E\2 

£*23 

Ux2 

18.844  (106) 

1.468  (106) 

0.91  (106) 

0.91  (106) 

0.45  (106) 

0.26 

Table  2.  Material  Properties  of  a  Hercules  AS4-3501-6  graphite  epoxy  com- 
posite(unit  psi) 


the  following  conditions;  conditionl:  dimensions  of  9  =1  radian,  24  plies 
[06/906]4  and  0.12  in.  thickness;  condition2:  9  =1  radian,  12  plies  [03/903]s 
and  0.6  in.  thickness;  ;  condition  3:  9  =0.5  radian  48  plies  [06/906]2s  and 
0.24  in.  thickness.  A  (176  nodes,  7x7  Elements)  quarter  symmetry  model 
was  used.  Results  are  shown  in  Figure  6,  7  and  8. 
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Load  (P/4)  (ibs) 


Displacement  (inch) 


Figure  6:  Cylindrical  Shell  (R=12in.  0  =  1  Radian,  24  Plies) 
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Figure  7:  Cylindrical  Shell  (R=12in.  6  =  0.5  Radian,  12  Plies) 


Figure  8:  Cylindrical  Shell  (R=12in.  6  —  0.5  Radian,  48  Plies) 
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From  Figure  6,  in  the  case  of  a  moderate  thickness  shell,  JAGS,  ABAQUS. 
and  LDLR  give  similar  results  up  to  and  beyond  collapse.  The  SLR  theory 
starts  to  depict  a  more  flexible  structure  in  the  unstable  region.  As  we 
described  in  theory  part  both  LDLR  and  SLR  are  Lagrangian  based  and 
employing  2nd  Piola-Kirchhoff  stresses  and  Green  strain.  The  difference  in 
SLR  theory  assumes  small  angle  approximations.  Because  the  thickness  isn’t 
large  enough  to  create  through  the  thickness  shear  effects,  the  resulting  value 
of  JAGS,  ABAQUS,  SLR  and  LDLR  are  quite  close.  From  Figure  7,  in  case 
of  the  very  thin  shallow  shell,  all  of  the  results  are  reasonably  close.  Finally, 
from  Figure  8,  for  the  shallow  and  thick  shell  case,  one  can  recognize  a  sig¬ 
nificant  difference  between  JAGS  and  the  other  programs  and  their  theories. 
This  difference  is  one  of  the  major  characteristics  of  the  JAGS  program.  It 
can  represent  through  the  thickness  shear  stresses  with  the  approximate  con¬ 
tinuities  better  than  any  of  the  other  theories.  SLR  doesn’t  include  a  through 
the  thickness  shear  stress  compatibility.  It  assumes  parabolic  shear  strain  ( 
it  assumes  third  order  shear  deformation).  LDLR  and  ABAQUS  assumes  1st  • 
order  shear  deformation  and  neither  of  these  theories  consider  direct  through 
the  thickness  shear  except  the  JAGS  theories  (based  on  corotational  theories, 
assumes  higher  order  shear  deformation). 


J 
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thickness 

E 

V 

radius 

0.0127  mm 

5.516  Gpa 

0.3 

1.5  m 

Table  3:  Geometrical  Properties  of  the  Shell  Structure  (Unit  psi) 

A  second  model  studied  was  a  circular  disk  (shown  in  Figure  9)  with  sim¬ 
ply  supported  boundaries  conditions  under  transverse  pressure.  This  prob¬ 
lem  has  been  solved  using  ABAQUS  and  compared  to  the  results  of  Greschik 
et  al.[4]  [5].  It  still  has  to  be  run  using  the  JAGS  program.  The  material 
properties  are  listed  in  Table  3. 


Finally,  the  third  model  considered  is  a  parabolic  shell  shown  as  a  dotted 
line  in  Figure  9.  (This  problem  still  has  to  be  rim  using  JAGS  and  was  only 
runs  using  ABAQUS.  The  results  were  compared  with  Greschik  et  al.  [4] 
[5]).  Figure  10  shows  930  S8R  type  elements  for  the  shell  model.  The  same 
material  properties  have  been  used  as  in  the  membrane  disk  except  it  follows 
parabolic  curves  ( it  depends  on  //£>;  where  /  represents  the  parabolic  focal 
length  and  D  is  diameter  of  parabolic  disk).  The  magnitude  of  zctr  charac¬ 
terize  the  shell’s  center  vertical  distance.  The  vertical  distance  zar  is  0  for 
flat  disk  and  7.376  inches  for  parabolic  membrane  shell  when  f/D  ratio  is 
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Figure  9:  Simply  supported  Disk  and  Parabolic  shell  under  transverse  pres¬ 
sure 

x2 

1.  The  equation  of  parabolic  disk  can  expressed  as  z  —  — -£■  —  zmt.  The 

(f) 

radial  and  circumferential  stress  components  for  node  93  (at  the  boundary 
near  edge)  under  an  increasing  internal  pressure  is  shown  in  Figure  11  for 
the  circular  membrane,  and  the  displacement  at  the  center  according  to  the 
increasing  pressure  is  showing  in  Figurel2.  One  should  note  the  nonlinear 
characteristics.  ABAQUS  does  have  the  capability  of  incrementally  increas¬ 
ing  the  pressure  and  maintaining  its  direction  normal  to  the  deforming  sur¬ 
face.  Thus  the  nonconservative  capability  was  implemented  in  running  the 
problems.  The  parabolic  shell  is  a  good  initial  depiction  of  the  antenna  sur¬ 
face,  and  thus  it  becomes  important  to  characterize  the  observed  results  for 
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low  loading  pressures.  Figures  13,  14  depict  the  stress  components  at  node 
93  for  the  parabolic  shell  geometry  and  the  center  displacement  in  which  the 
supports  were  assumed  to  be  pinned  (radial,  tangential  and  z  direction  is 
constrained  as  fixed).  These  two  figures  show  that  the  support  condition 
creates  the  type  of  stress  field  recognized  for  its  linear  membrane  charac¬ 
ter.  The  tangential  stress  is  linearly  increasing  and  is  continuously  showing 
a  tension.  The  center  deflection  is  also  purely  linear  relative  to  the  load  as 
shown  in  Figure  14.  In  the  next  condition  considered,  the  radial  constraint 
is  removed.  The  deep  shell  with  a  / /D  =  1  indicates  that  a  compressive 
stress  component  occurs  in  the  tangential  or  hoop  direction  and  Figures 
15  -  18  are  verification  of  this  stress.  At  node  185  negative  hoop  stresses 
increase  when  the  loading  increases.  The  location  of  nodes  93  and  185  are 
described  in  the  Figure  10.  Figures  19  -  22  represent  the  fact  that  radial  and 
hoop  stresses  change  according  to  coordinate  angle,  (where  the  coordinate 
angle  is  measured  counter  clockwise  with  respect  to  the  x-axis;  x-axis  is  0 
and  y-axis  is  90  degree).  Negative  hoop  stresses  exist.  It  can  also  be  noticed 
that  close  to  the  supports  the  compressive  hoop  stress  region  decreases  with 
respect  to  a  load  increase  as  shown  in  Figures  21  and  22.  In  fact,  at  a  load 
of  5.9531  x  10  4  psi,  the  shell  shows  a  very  small  value  of  compressive  hoop 
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stress.  Figure  23  and  24  represent  the  hoop  stress  value  along  a  coordinate 
angle  of  close  to  zero  degrees.  One  can  observe  the  region  of  negative  val¬ 
ues  of  this  stress  component  and  their  magnitude  at  load  1.2711  x  10-4  and 
3.1777  x  10-4  psi.  It  becomes  apparent  that  the  boundary  restraint  is  very 
important  in  relationship  to  the  membrane  stress  component. 
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coordinate 

angle 

measurements 


Figure  10:  Node  numbers  and  Element  numbers  of  Flat  disk  and  Parabolic 
membrane  shell  geometry,  (regular  numbers  represent  node  numbers  and 
underlined  numbers  represent  element  numbers)  Representation  of  Radial 
stress  Srr  and  Tangential  stress  Stt. 
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Stress(psi) 


Figure  11:  Radial  Stress(Srr)  and  Tangential  Stress(Stt)  vs.  Load  Curve 
( Initially  flat  disk  with  pinned  support  case). 
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Load  (lbs) 


Figure  12:  Load-Displacement  Curve  ©Center  (Initially  flat  disk  with 
pinned  support  case). 
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Load  (lbs) 


Figure  13:  Radial  Stress(Srr)  and  Tangential  Stress(Stt)  vs.  Load  Curve 
(  Initially  Parabolic  membrane  shell  geometry  (f/D  =  1)  with  pinned  sup¬ 
port  case) 
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Load  (lbs) 


Figure  14:  Load-Displacement  Curve  ©Center  (Initially  Parabolic  mem¬ 
brane  shell  geometry  (//£>  =  !)  with  pinned  support  case). 
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Srr,  Stt  vs. Load  Curve  @node=93 


Figure  15:  Radial  Stress(Srr)  and  Tangential  Stress(Stt)  vs.  Load  Curve  @ 
node  93  ( Initially  Parabolic  membrane  shell  geometry  (f/D  =  1)  with 

no  radial  constraint-pinned  support  case). 
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Stt  vs  Load  @node=93 


Figure  16:  Detailed  figure  of  Tangential  Stress(Stt)  vs.  Load  Curve  @  node 
93 

(  Initially  Parabolic  membrane  shell  geometry  ( f/D  —  1)  with  no  radial 
constraint-pinned  support  case). 
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Sit,  Stt  vs.  Load  Curve@node=185 


Figure  17:  Radial  Stress(Srr)  and  Tangential  Stress(Stt)  vs.  Load  Curve  @ 
nodel85 

(  Initially  Parabolic  membrane  shell  geometry  ( f/D  =  1)  with  no  radial 
constraint-pinned  support  case). 
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Stt  vs.  Load  @node=185 
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Figure  18:  Detailed  figure  of  Tangential  Stress(Stt)  vs.  Load  Curve 
(  Initially  Parabolic  membrane  shell  geometry  (f/D  =  1)  with  no  radial 
constraint-pinned  support  case) 
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Srr,  Stt  vs.  Coordinated  Angle  @r«47.581  inch 


Figure  19:  Radial  Stress(Srr)  and  Tangential  Stress(Stt)  @  radial  dis- 
tance=47.581  inch  and  load=  1.2711  x  10-4  psi 

(  Initially  Parabolic  membrane  shell  geometry  (f/D  =  1)  with  no  radial 
constraint-pinned  support  case). 
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Negative  Stress  Location 

Srr,  Stt  vs.  Coordinated  Angle  @r=49.484  inch 


-20 


Coordinate  Angle  (degree) 

Figure  20:  Radial  Stress(Srr)  and  Tangential  Stress(Stt)  @  radial 
distance=49.484  inch  and  load=  1.2711  x  10-4  psi 

(  Initially  Parabolic  membrane  shell  geometry  (f/D  =  1)  with  no  radial 
constraint-pinned  support  case). 


3° 

Negative  Stress  Location 


Srr,  Stt  vs.  Coordinate  Angle  @r=49.484  Inch 


Figure  21:  Radial  Stress(Srr)  and  Tangential  Stress(Stt)  @  radial  dis- 
tance=49.484  inch  and  load=  3.1777  x  10“4  psi 

(  Initially  Parabolic  membrane  shell  geometry  (f/D  =  1)  with  no  radial 
constraint-pinned  support  case). 
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1° 

Negative  Stress  Location 


Srr,  Stt  vs.  Coordinate  Angle  @  r=49.484  Inch 


Figure  22:  Radial  Stress(Srr)  and  Tangential  Stress(Stt)  vs.  Loads  Curve  @ 
radial  distance=49.484  inch  and  load=  5.9531  x  10-4  psi 
(  Initially  Parabolic  membrane  shell  geometry  (f/D  =  1)  with  no  radial 
constraint-pinned  support  case). 
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Stt  vs.  radial  direction 
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Figure  23:  Tangential  Stress  vs.  Radial  Distance  @  load  =  1.2711  x  10-4  psi. 
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Stt  vs.  radial  distance  @Load=3.1777x104  lbs 


Figure  24:  Tangential  Stress  vs.  Radial  Distance  ©  load  =  3.1777  10 ~4psi. 
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4  Further  Research 


This  effort  is  directed  toward  the  development  of  a  program  that  is  specifically 
designed  to  carry  out  the  nonlinear  analysis  of  antenna  structures.  JAGS  is 
the  in-house  program  worked  upon  for  this  purpose.  So  far  we  have  char¬ 
acterized  an  individual  panel,  a  disk  and  a  parabolic  shell  structure  under 
various  loadings  using  ABAQUS.  The  JAGS  program  has  been  compared  to 
ABAQUS  for  the  panel.  The  equations  have  been  developed  for  nonconserva¬ 
tive  loading  to  be  implemented  in  JAGS.  Also,  the  transformation  relations 
needed  for  a  generalized  global  set  of  elements  was  developed  and  needs  to  be 
implemented  in  JAGS.  Once  the  new  relations  for  a  nonconservative  loading 
and  a  generalized  transformation  are  implemented  in  JAGS,  the  membrane 
and  parabolic  shell  will  be  run.  The  authors  will  also  continue  to  characterize 
the  antenna  structure  using  ABAQUS.  We  will  next  combine  two  parabolic 
shells;  one  concave  up  and  the  other  concave  down  joined  at  their  connecting 
points  under  internal  pressure.  This  will  be  followed  by  including  discrete  tie 
rods  at  specific  locations.  The  final  model,  that  is  the  global  configuration, 
will  be  one  which  includes  the  parabolic  shell  configurations  joined  to  a  torus 
through  the  appropriate  tie  rod  configurations.  The  torus  will  be  supported 
a  three  concentrated  points  separated  by  120  degree.  It  is  apparent  that 
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the  most  important  part  of  the  last  scenario  is  the  instability  of  the  torus. 
The  JAGS  program  results  will  be  compared  to  ABAQUS  throughout  the 
research.  As  the  global  model  becomes  more  and  more  complicated,  it  has 
become  obvious  that  the  computer  intensive  feature  will  be  a  challenge. 
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Hi  =  (1/8)  (1  4-  rkr)  (1  +  sks)  (2  +  rkr  +  sks  -  r2  -  s2)  (A3a) 

H\  =  (a/8)  rk  (1  +  rkr)2  (rkr  —  1)  (1  +  sks)  (A3b) 

Hkz  =  (b/8)rk(l  +  rkr)(sks-l)(l  +  sks)2  (A3c) 
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Ck  =  (1/ 4)  (1  +  Tkr)  (1  +  SkS) 


(A3d) 


where  2a  and  2b  are  dimensions  along  x  and  y  of  the  rectangular  (in  curvi¬ 
linear  coordinates)  element,  and  the  values  of  and  s*.  are  determined  by 
the  local  coordinates  (r,s)  of  the  kth  node. 

A  vector-^}  displacement  quantities  at  the  reference  surface  can  be  ex¬ 
pressed  as; 

{ip}  =  {(1  +  ei)cos761  —  1,  (1  +  e2)  cos762,  (1  +  e!)sin761 -f  (1 -(- e2)  sin762, 

ki  -  k2  -  k%,  ke  -  k%,  74)!E,  74  y,  /y5x,  75}T  (A4) 

{ipo}  3X6  all  known  values  at  q0  in  equation(18). 

A  expression  12x12  symmetric  matrix  [$>]  is 

m  =  EjT”  [S“»]T  3W]  [sw]  <fe  (A5) 

where  Qy  is  the  material  transformed  stiffness  matrix  for  the  ith  lamina 
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from  its  principal  matrix  and  the  matrix  is 
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(A6a) 

The  values  of  the  entries  g$  (warping  and  thickness  stretching  function)  are 

#51 

)  _ 

fagS  +  k5942, 

952  —  ^59 14  +  fag24 

(A6aa) 

#61  — 

h943  +  fa9u> 

962  —  fa9l5  +  ^4^25 

(A6ab) 

9% 

)  = 

a3j  +  2^4- 

3 

=  0,..,9; 

9«  =  a«  +  2$ 

j 

=  1,  ..,4  (A6ac) 

9$ 


n (0 
954 


924, z  ~  {fa924  +  fa29u)  ,  963  ~  925, z  ~  (fa925  +  ^625l5)(A6ad) 
9u,z  ~  (fai924  +  kiffu) ,  g$l  =  <7j5]2  —  (^6x^25  +  fciyis)(A6ae) 


An  dement  of  the  [4/]  matrix 

To  illustrate  the  content  of  the  [\J>]  matrix,  we  first  consider  the  displace- 
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ment  gradient  vector,  {U} 


{t/}  {u,  u,x  ;  u,y ;  v ;  v,x  ;  v,y  ;  v)XX  ;  viXy  ;  v,yy ; 


24x1 


W,W1X  ,W,y  ,  W^XX  ,  t*L) jXy  ,  5  74!  74  5  74?!/  5 

75:7s*  !  }T 


(AT) 


and  the  variation  of  the  vector 


<5  {^}  =  [^]<5{C/} 


(AS) 


where 


%  =  [¥]  = 


8% 

dUj 


(A8a) 


’*r  (iJ)  functions  of  the  global  displacements  u,  v,  w  and  the  initial  curva¬ 
ture.  The  terms  are  generated  in  mathematica.  It  can  be  shown  for  example 
d% 


that  ¥(6,4)  = 


du,xx 

¥(6,4)  =  dip6/dUA 

=  d  (k6  -  /du,xx 
—  9  (k$\  +  k$2  —  —  k fdu 


)XX 


(A9) 


Q»  O* 

where  A;61and  fc62  are  the  scalar  quantities  representing - -  •  i2  and  -  —  •  i3 

dx  dy 

respectively  (see  Figure  2),  and  A&  for  example  is  the  curvature  term  that 
yields  the  j3  component  of  the  again  see  Figure  2). 
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The  results  is 


®  (6)  4)  —  d  (— T2\yXTz\  —  T22,xTs2  ~  T22,xT$z)  / du,xx  (A10) 

where  the  axe  functions  of  fiber  stretching  and  the  corotational  axes  at 
the  deformed  configuration  of  a  point.  This  yields  the  expression  for  \fr  (6, 4) 
in  terms  of  a  numerator  and  denominator  function; 

*(6,4)„  =  c  (—  (ab)  +  i>2  +  d2) 

x •fa  (a5  -  2a36  +  aS2  +  Ztfd*  -  bd?  -  MP  -  2\/a<P\/5^'35) 

X  (-  (C3C5)  +  C2Ce)  (AlOa) 

^(®>4)d  =  \/a\/ (°3  —  ab  +  d2)2  (b  +  y/ab  —  d2) 

xjb  (a26  -  2ab2  +  ft3  -  ad2  +  3bd?  -  adVafc  -  d2)  (AlOb) 

X  +  C3^  ~  2C1C2C4C5  +  C1C5  +  C3C5  —  2C1C3C4C6  +  cfcg  +  C^Cjj 

where 


Ci  =  1  +  u,x  —vk®  + 1 ok° 

C2  =  v,x  — ufc®  +  tcfcg! 

C3  wfcj  “b 

C4  —  Ujy  ~~vk^  xvkfi 2 
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cs  -  1  +  v,y  -uk°  +  Wk% 
c«  =  u/.„  -U^2  +  vk$ 


(AlOc) 


a  ~  yfk  +  C2+<^ 

^  =  v^Ti 

C  =  C1C4  +C2C5  +C3C6 

d  =  i/v^6 

All  other  terms  are  expressed  in  detail  in  Appendix  of  Creer[14]  or  refer 
to  Greer  and  Palazotto[15], 
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